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In this paper we present results from numerical calculations for three identical boson systems for 
both very large and infinite two-body s-wave scattering length a. We have considered scattering 
lengths up to 2 x 10 J a.u. and solved the hyperangular part of the Schrodinger equation for distances 
up to I0 6 a.u.. Form these, we obtained the three-body effective potentials, hyperspherical channel 
functions and the asymptotic behavior of the nonadiabatic couplings in order to to characterize the 
main aspects of the Efimov states. These results allow us to test and quantify the assumptions 
related to the Efimov effect. 
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I. INTRODUCTION 



The experimental ability to create ultracold atomic 
systems, such as Bose-Einstein condensates, with a large 
range of interaction strengths has raised the possibility 
of investigating a very unusual phenomenon: the Efimov 
effect. The Efimov effect |l| is responsible for the emer- 
gence of an attractive long-range effective three-body in- 
teraction irrespective of the short-range character of the 
atomic interactions, producing a large number of three- 
body bound states — called Efimov states — even when 
the two-body subsystems have none. This effect is ex- 
pected to occur in the limit where the two-body s-wave 
scattering length a is much larger than the characteristic 
range of the interatomic interactions. Such a regime 
can be reached, for instance, by tuning the two-body in- 
teractions via an external magnetic field near a Feshbach 
resonance Li LjJ- Li LJl ■ Besides its intriguing impact 
on the three-body bound state spectrum, Efimov physics 
has been invoked to explain the remarkable universality 
of three-body scatte ring observables at ultracold ener- 
gies 0EUH[iJ[ilH0. Here > we use the less 

restrictive term "Efimov physics" to indicate the quali- 
tative change in behavior exhibited by the system when 
\a\ — > oo. 

Despite their importance, neither Efimov states nor 
any other signature of Efimov physics has yet been de- 
tected. The main difficulty for direct observation of 
these states stems from the fact that the Efimov states 
are very weakly bound and easily perturbed. Indirect 
observation of Efimov physics might be possible, how- 
ever, by measuring some three-body scattering observ- 
ables like three-body recombination or vibrational relax- 
ation 0, H, El El El as a function of the scatter- 
ing length. The detection of the series of interference 
minima and resonance peaks present in the three-body 
recombination rate as a function a, or of the resonance 
peaks in the vibrational relaxation rate, would be a clear 
signature of the Efimov physics involved in the ultra- 
cold three-body collisions. In this case, the main diffi- 
culties arise from the fact that finite collision energies 
- even ultracold energies — reduce the range of a in 
which Efimov physics can be observed [T3. At finite 



energies, the unitary limit, high partial waves, and ther- 
mal averaging tend to wash out the features connected to 
Efimov physics. The coldest current experiments, how- 
ever, might offer a large enough window in scattering 
length to see Efimov physics, but they would still be dif- 
ficult experiments. For Bose-Eistein condensates of 4 Hc, 
23 Na and 87 Rb, for instance, we can expect to observe 
at most three interference minima or resonance peaks in 
the three-body recombination rate for energies of 3 nK, 
0.5 nK and 0.1 nK, respectively To see more than 
three of either requires even lower temperatures. 

Previous theoretical work 0, H H 0, IH E2 has re- 
vealed many universal low-energy three-body scattering 
properties related to Efimov physics. In this regime, the 
effective interactions and the three-body observables are 
determined primarily by a, independent of the details of 
the two-body interactions. The assumption that the at- 
tractive long-range behavior of the three-body effective 
potential (ex —1/R 2 ) is valid in the range r -C R <C |a|, 
with R giving the overall size of the three-body system, 
is usually required in order to determine the a depen- 
dence in the three-body obsevables. In fact, the more 
relaxed assumption that this potential holds even for 
tq ^ R H has been used to deduce the number of Efi- 
mov states and other low-energy three-body observables 
P. Hol ITU IT^|. Neither assumption, however, has been 
thoroughly tested against accurate calculations. Given 
the central role these assumptions play in determining 
the universal properties of ultracold three-body scatter- 
ing, they should be tested. 

For this paper, we have solved the Schrodinger equa- 
tion for three identical bosons using the adiabatic hyper- 
spherical representation for the J n = + symmetry over 
a broad range of scattering lengths, including essentially 
infinite values. The numerical results were obtained in 
a regime not accessible in previous calculations, allowing 
us to test the assumptions above by checking the modifi- 
cations to the effective three-body potentials as \a\ — > oo. 
In addition, we have explored the validity of the above 
assumptions as the number of two-body bound states in 
the two-body potential is increased. We have consid- 
ered both positive and negative scattering lengths, up to 
2 x 10 5 a.u. in magnitude and separately the infinite a 
case. We solved the three-body Schrodinger equation for 
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distances R as large as 10 6 a.u. in order to determine 
the effective three-body potentials as well as the nonadi- 
abatic couplings responsible for inelastic transitions. The 
asymptotic behavior of the coupling and the hyperspher- 
ical channel functions was investigated to point out the 
main characteristics of the Efimov states in the adiabatic 
hyperspherical point of view. 



II. ADIABATIC HYPERSPHERICAL 
REPRESENTATION 

We solve the three-body Schrodinger equation using 
the adiabatic hyperspherical representation 0, 0] . Af- 
ter the usual separation of the center-of-mass motion, the 
three-body problem can be described by the hyperradius 
R and five angles, denoted collectively by fl. The five an- 
gular coordinates are chosen to be the Euler angles (a, 
(3 and 7) that specify the orientation of the the plane 
defined by the three particles relative to the space-fixed 
frame, plus two hyperangles (ip and 9) defined as a mod- 
ification of Smith- Whitten coordinates [TEl HH . Il7| . 

The Schrodinger equation in hyperspherical coordi- 
nates can be written in terms of the rescaled wave func- 
tion tp = i? 5 / 2 ^/ (in atomic units), 



1 d 2 

^'dR 2 



H, d (R,n) 



^(R,n) = Ei,(R,n), (1) 



where fi is the three-body reduced mass (we choose /j, 



m/y/3 for three identical bosons of mass m) and E is the 
total three-body energy. The adiabatic Hamiltonian H ad 
is given by 



H ad (R,n) 



A 2 (») + f 
2fiR 2 



(2) 



in which A 2 is the grand angular momentum operator 
and V contains the interparticle interactions. 

In the adiabatic hyperspherical representation, the to- 
tal wave function is expanded on the channel functions 



ii(R,n)=J2F„(R)$ 1/ (R;n), 



(3) 



where the expansion coefficient F V {R) is the hyperra- 
dial wave function. The total wave function is still, in 
principle, exact. The channel functions & V (R\ f2) form a 
complete set of orthonormal functions at each R and are 
eigenfunctions of the adiabatic Hamiltonian: 



H ad (R, fl)<f> v (R; 0) = U V (R)$ V {R; fi). 



(4) 



The eigenvalues of this equation, U V {R), play the role 
of effective potentials for the hyperradial motion. This 
fact can be seen upon substitution of Eq. J3J into the 
Schrodinger equation JQ| and projecting out $„. The re- 
sult is a system of coupled ordinary differential equations 



1 d 2 
2n dR 2 



U V {R) 



F„{R) 



-E 



2P VV \R) 



dR 



Qvu'{R) 



F V ,(R)=EF V (R), 



(5) 



where P UV '(R) and Q UV '(R) are the nonadiabatic cou- 
pling terms. These couplings are generated by the in- 
dependence of the channel functions and are responsi- 
ble for the inelastic transitions in three-body scattering. 
They are defined as 



and 



®u(R) 



dR 2 



(6) 



(7) 



where the double brackets denote integration over the an- 
gular coordinates only. Here, we will concentrate the 
solutions of the adiabatic equation @, namely the poten- 
tial curves and couplings. For the present calculations, 
the potential V in Eq. J2J) is a pairwise sum of short-range 
potentials, V — u(ria) + v(r3%) + v(r23), which is a good 
model for spin-stretched atoms. 



A. Two-body potential model 

To a good approximation, the Efimov effect, as well as 
other low-energy three-body properties, do not depend 
on the details of the interatomic interaction, only on the 
scattering length and ro, the characteristic range of the 
interatomic interactions. This approximation, verified for 
different two-body potentials 0] , allows the use of model 
potentials, substantially simplifying the numerical calcu- 
lations by controlling the number of two-body channels 
in the problem. The two-body potential adopted here is 
v(r) = Dsech 2 (r / tq) . With this potential, we can easily 
produce two-body systems with different values of the 
scattering length by changing the interaction strength 
D. The scattering length is determined for each value 
of D by solving the two-body scattering problem. For 
infinite values of the scattering length, however, we used 
a slightly different approach which will be discussed in 
more detail below. We have considered potentials with 



3 



just one two-body bound state and a characteristic range 
of r = 15 a.u., and using the mass of 4 He atoms. 

Figure shows the scattering length a as a function 
of the potential strength D. Each time that the two- 
body potential deepens enough to support an additional 
s-wave bound state, the scattering length goes through a 
pole, labeled in Fig.^by Dj, Du and Dm. For D > Di, 
the two-body system has no bound state and represents 
the standard Efimov case when \a\ — > oo. For Du < 
D < Dj, the system has a single s-wave bound state, and 
for Djjj < D < Djj, the system has gained another s- 
wave bound state. Our model two-body potential mimics 
the magnetic field control of the scattering length near 
a two-body s-wave Feshbach resonance Q, 0, m IE 111 
assuming that the only effect the resonance is to change 
the scattering length. In fact, we expect this model to 
be valid in the threshold regime, i.e., when the collision 
energy is the smallest energy in the system As a -> 
+oo, the two-body bound state is weakly bound, and the 
three-body effective potential that approaches this state 
asymptotically supports the three-body Efimov states. 
As a — > — oo, the two-body bound state becomes a deep 
bound state and the Efimov states are resonances rather 
than bound states. 
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FIG. 1: The two-body s-wave scattering length a as a function 
of the potential strength D. When the two-body potential is 
deep enough to support an additional bound state, a passes 
through a pole, indicated by Di, Du and Dm. 

To determine the value of D that gives a particular 
scattering length, we simply read the value from Fig. ^ 
To determine the value of D that gives \a\ — > oo from the 
figure, though, is not simple to do accurately. We thus 
used a slightly different approach that takes advantage of 
the fact that the derivative of the two-body radial wave 
function must vanish at large distances when \a\ = oo. 
To this end, we solved the two-body radial Schrodinger 
equation at zero energy, 



1 d 2 

— -r~9 + Dsech 2 (r/r ) 
m dr z 



u(r) = 0, 



(8) 



and required u'(r) = at some large distance. The po- 
tential strength D is the eigenvalue of Eq. |JSJ — the low- 



est corresponds to the solution with no nodes D = Dj, 
the next to the solution with one node D — Du, and so 
on. This procedure determines D to essentially machine 
precision. We can then perform three-body calculations 
with a scattering length as close as numerically possible 
to infinity. 

It is worth to noting that the solutions of the sech 2 
potential are known analytically, so that we could, in 
principle, determine these values of D exactly. We used 
the present method because it will work for any short- 
range potential. 



B. Effective hyperspherical potentials 

Neglecting all off-diagonal terms in Eq. (J5J defines the 
effective hyperspherical potential, 



W U {R) = U V {R) - 2jfi*u(R). 



(9) 



For short-range two-body potentials the channel index v 
labels both two-body and three-body configurations in 
the asymptotic limit (R — > oo). In this limit, the po- 
tential curves U V (R) approach either a two-body bound 
state energy for the bound channels (BC) or zero energy 
for the three-body continuum channels (CC). In both 
cases, Q vv vanishes for R — > oo. 

For finite a, however, the effective potentials arc influ- 
enced by Efimov physics in the range r^<^.R<^]a\. For 
these R, the effective potentials areproportional to R~ 2 , 
and can be attractive or repulsive. The channel of prim- 
iry interest in this case is the one that leads to the Efimov 
effect — for a > 0, it is the chanel correlated asymptoti- 
cally with the weakly bound s-wave dimer; for a < 0, it 
is the lowest continuum channel. Because of its impor- 
tance in our discussion we label it the "Efimov channel" 
(EC). It is worth to mention that for R ^> \a\ the Efi- 
mov channel behaves as a bound or continuum channel, 
as appropriate. 

In Fig. |21 we show the effective potentials W V (R) cal- 
culated near the pole D = Dj in Fig. ^ For a > 0, the 
lowest potential is the Efimov channel (EC), converging 
to the two-body bound state for R — > oo. The higher ly- 
ing potentials are associated with three-body continuum 
channels (CC). As D is increased towards Dj from the 
left, a — ► +oo and the dimer binding energy vanishes. 
At the same time, the Efimov channel evolves from that 
shown in Fig. [21 to a purely attractive R~ 2 for R ^> r*o- 
Similarly, for a < in Fig. |2{b), the lowest potential is 
the Efimov channel (EC) converging to the three-body 
continuum limit as R — > oo. In this case, a — ► — oo as 
D approaches Di from the right, and the potential bar- 
rier in the Efimov channel moves to infinity (R oc \a\). 
At D = Di, the potential is the same attractive R~ 2 
obtained by approching Dj from the left. Note that for 
both a < and a > 0, there can be bound channels (BC) 
that lie at lower energies. 
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FIG. 2: The three-body effective potential for (a) a > and 
(b) a < 0, calculated near pole D — Di. For a > 0, the 
system has only one bound state, while for a < the system 
has none. 



The asymptotic behavior, R 3> \a\, of the effective 
potential W v (R) can be determined analytically 0, 0] . 
For the bound channels (BC), the asymptotic limit repre- 
sents the atom-dimer configuration, and the asymptotic 
effective potential is given by 



W V {R) = E vl 



1(1- 



2fiR 2 



(10) 



The rovibrational energy E v y of the dimer is labeled by 
its vibrational quantum number v and its angular mo- 
mentum V; I is the relative angular momentum between 
the dimer and the third atom. The notation C(...) indi- 
cates the leading order correction. 

The three-body continuum channels (CC) represent 
three free atoms, each one far from the others in the 
asymptotic limit. The two-body potentials thus vanish, 
leaving only the kinetic energy. Therefore, the leading 
term of the effective potential is determined by the eigen- 
values A(A + 4) of the kinetic energy, and the higher order 
terms can be written in terms of the scattering length a, 



W V {R) = 



A(A + 4) 



15 
4 



2nR 2 
C2 f a\ 2l «+ 2 



[iR 2 \R 



liR 2 

+ o 



a \ 2*0+1 



(I) 



1 



where c\ and c-i are positive constants. We can relate the 
quantum number A to the partial angular momenta of the 



system using A = 2n + l + l' — n is a non-negative integer, 
V is the angular momentum associated with the Jacobi 
vector connecting particles 1 and 2 and I with the Jacobi 
vector connecting particle 3 with the center of mass of 1 
and 2 [see also Eq. (|10l) ] . For identical bosons with J 71 — 
+ , symmetry requires A = 0, 4, 6, and I = V = 0, 2, 
4, .... The quantum number Iq in Eq. (| 1 1 ft can thus be 
thought of as the dominant partial wave in the expansion 
of a continuum state labeled by A. It is non-zero only 
when there are degeneracies in A — which begin in the 
present case with A = 12. Thus, this situation does not 
occur for the lowest continuum channel which is usually 
of most interest for ultracold collisions. We will use A 
and Iq to label the three-body continuum channels 

As discussed above for \a\ — > oo, the effective potentials 
for both a < and a > in the range ro<^R<^\a\ are 
proportional to R~ 2 . Explicitly, 



W V {R) 



2[iR 2 



O 



1 

R* 



(12) 



with so — 1-00624, therefore producing an infinite num- 
ber of three-body bound states when \a\ = oo — which is 
precisely the Efimov effect. The higher- lying continuum 
channels also feel the influence of the Efimov physics and 
are no longer associated simply with the eigenvalues of 
the kinetic energy. For ro<Ci?<C|a|, they are given by 



W V (R) 



2[iR 2 



O 



1 



(13) 



where s„ is a constant that represents a sort of correlation 
in addition to the kinetic energy. 

In some way, the assumption that the potentials l|12l) 
and hold for r < R < \a\, has been used [9U10 UIII 
llH to deduce many properties of three-body systems with 
a large scattering length. In fact, it was the more relaxed 
assumption that the long-range effective potential l|12l) 
holds even for ro < R < \a\ that was used in Refs. pi ITd 
0,0] to deduce some low-energy three-body observables 
and the number of Efimov states, 



N ~ (s /n)ln(\a\/r ). 



(14) 



The validity of the assumption that ro < R < \a\ sets 
limits on the quantitative aspects of the these predic- 
tions. In the next section, the above assumptions will be 
tested with accurate numerical calculations for scattering 
lengths up to 200000 a.u. (and tq — 15 a.u.), which is 
unambiguously in the regime \a\ Tq. 



III. RESULTS AND DISCUSSION 



The adiabatic equation 10} has been solved for hyper- 
(11) radii up to 2 x 10 6 a.u. in order to reach the asymp- 
totic regime for the effective potential and couplings. We 
can thus compare our numerical results with the com- 
mon assumptions of the Efimov effect. The main thrust 
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of this investigation is to establish quantitative bounds 
on the range where the Efimov potential accurately ap- 
proximates the effective potential. In addition, we have 
investigated the qualitative differences between Efimov 
states and ordinary three-body bound states by analyzing 
the channel function and nonadiabatic couplings when 
\a\ — > oo. 



A. Effective Potentials 

We show in Fig. the effective potentials associated 
with the weakly bound dimcr for a > 0, and with low- 
est continuum channel for a < 0, for several values of a. 
They are multiplied by the factor 2/ii? 2 in order to more 
clearly reveal whether they match Eq. Ijl2|l . If Eq. HI2(I 
holds, then these curves should approach the constant 
— Sg(= 1.0125). Figure shows our numerical potentials 
for \a\ — oo (solid lines) and finite a (dashed lines) as 
well. For a > in Fig[3Ja), the effective potentials cor- 
respond to the weakly bound channel, converging to the 
energy of the weakly bound dimer state [Eq. (|I0|l ] given 
approximately by — I /ma 2 . This fact manifests itself at 
large R in Fig. [31 by the — R 2 divergence of the curves. 
As the scattering length grows in magnitude, the effec- 
tive potential should converge to the Efimov potential 
(indicated by the horizontal dashed line) in the range 
r < R < |o|. 

Similarly, Fig. UJb) shows the a < effective poten- 
tials corresponding to the lowest continuum channel and 
converging to the Efimov potential when \a\ — > oo (hori- 
zontal dashed line). For finite a, this curve must asymp- 
totically approach the lowest kinetic energy eigenvalue 
(dotted line) — A(A + 4) + 4 with A = 0. The lowest 
effective potential converges to the deeply bound dimer 
state and is omitted. Efimov states in this case are thus 
resonances rather than bound states, since they can de- 
cay to the lower-energy dimer channels Q . 

The vertical lines in Fig. UJ correspond to R — r$ and 
R = |o|, with \a\ = 500, 10000, 100000 and 200000 a.u.. 
Within this region — actually ro < i? < \a\ — the 
associated potential is expected to be universal, i.e., be- 
have like the Efimov potential. Clearly, the potentials 
do not satisfy this expectation over the whole range of 
ro < R < | a | as is often assumed. Instead, the condition 
specifies R "much" greater than ro and "much" less than 

14 

If we arbitrarily establish a limiting relative error of 
10% in the comparison between the finite scattering 
length potential, as plotted in Fig. [3J and the Efimov po- 
tential, then we see that for positive scattering lengths in 
the range 20ro < R < |a,| / 13 the effective potentials can 
be considered an Efimov potential. For a < 30000 a.u., 
however, the effective potential does not come within 10% 
of the Efimov potential in the range ro < R < a. For 
negative scattering lengths, the effective potentials can 
be considered (with an error of 10%) to be an Efimov 
potential in the range 33ro < R < |a|/13. Thus, both 
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FIG. 3: Three-body effective potentials for (a) a > and 
(b) a < 0, with |o| = 500, 10000, 100000, and 200000 a.u. 
(dashed lines). Also shown is the numerical result for \a\ — oo 
(solid line). The horizontal dashed line at — Sq is the Efimov 
potential, Eq. (11 21 , and the vertical lines correspond to R — ro 
and R = \a\. 



positive and negative a results are reasonably consistent 
with the assumption that the Efimov result holds for 
r <C R <§; \a\, but are not consistent with r < R < \a\. 
Since the standard estimate for the number of Efimov 
states, Eq. I|14(l . assumes that the Efimov potential holds 
in the full range of R from ro to \a\, these more restrictive 
ranges reduce the estimated number of Efimov states by 
about 2. 

For the 4 He three-body problem, for instance, where 
a ~ 170 a.u., the effective potential does not approach 
within 10% of the Efimov potential in the range ro < R < 
a. The above-mentioned restrictions lead to the conclu- 
sion that no Efimov state should occur in this case. Nev- 
ertheless, accurate numerical calculations predict a single 
Efimov state [2(j, which agrees with the result obtained 
using Eq. I|14f> where no such restrictions were imposed. 
Therefore, despite the more restrictive assumptions dis- 
cussed here, Efimov physics can still be observed — even 
when the effective potential does not approach the Efi- 
mov potential. The potentials are still modified in the 
range r < R < \a\, just not as dramatically. This modi- 
fication can be seen in Fig.|21for a — 500, where the curve 
shows a pronounced maximum near R = 100 a.u. before 
assuming its asymptotic behavior. This sort of deviation 
seems to indicate the influence of Efimov physics — the 
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smaller this feature, the weaker the influence. 

In Fig. 0] we show the effective potential for the Efimov 
channel with \a\ — oo for systems with different num- 
bers of s-wave bound states. The labels Dj, Du, Dm 
and Djv denote the two-body potential strengths that 
produce \a\ = oo with no s-wave bound states, and with 
one, two and three s-wave bound states, respectively (see 
Fig. ^) . For systems with multiple s-wave bound states 
there are also rotationally excited two-body bound states 
that are included in the numerical calculations. The inset 
in Fig.0|shows that the lower limit on the range of valid- 
ity of the Efimov potential (|f 2fl changes with the number 
of two-body bound states. Therefore, the restriction on 
the lower limit grows even stronger when the number of 
bound states increases. For instance, with two s-wave 
bound states {Dm in Fig. the lower limit must be 
set at 48ro to meet our 10% relative error criterion. Fig- 
ure 01 also appear to converge to a limiting curve as the 
number of two-body bound states increases. 
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FIG. 4: The attractive dipole effective potential for different 
number of two-body bound states. Di, Du, Dm (as given in 
Fig. and Div (not shown in Fig. indicates the effective 
potentials with no s-wave bound states, and with one, two 
and three s-wave bound states, respectively. 

It is interesting to notice that for \a\ — > oo the con- 
tinuum channels also reflect Efimov physics. While they 
still represent three free particles, they also include some 
sort of correlation in addition to the kinetic energy. Fig- 
ure shows the effective potentials for the continuum 
channels with \a\ — oo (solid lines) and a = 8000 a.u. 
(dashed lines). As R — > oo, the curves for finite a ap- 
proach the eigenvalues of the kinetic energy A(A + 4) + 4 
(the dotted lines in Fig. [5J). The effective potentials for 
\a\ = oo, however, converge to a different limit when 
R — > oo, and can be determined analitically pj. From 
Fig. we can see that some of the curves for \a\ = oo 
still converge to the kinetic energy eigenvalues. Gener- 
ally, these curves are degenerate in A with l ^ — only 
states that have primarily s-wave two-body character are 
impacted by Efimov physics. As was the case for the at- 
tractive dipole potential, the effective potentials for the 
continuum channels with finite a behave like the Efimov 



potentials [Eq. I|f 3(1 ] in the range ro <C R <C \a\. In fact, 
these modifications to the effective potentials have been 
seen for other three-body systems (also for higher partial 
waves) and have been shown to play a fundamental role 
in the scattering length scaling laws of ultracold three- 
body collisions [l3| . 
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FIG. 5: Three-body effective potentials associated with the 
continuum channels for \a\ = oo (solid lines) and for a — 
8000 a.u. (dashed lines). The dotted lines correspond to the 
eigenvalues of the kinetic energy; and the vertical dot-dashed 
lines to the boundaries ro and a. 



B. Channel Functions 

In Fig. we show the channel functions for the three 
classes of states: bound (BC), continuum (CC), and Efi- 
mov (EC). For the J n = + symmetry, the channel func- 
tions do not depend on the Euler angles; they can thus be 
visualized for fixed values of R as a function of the hyper- 
angles (f and 8 only. We show only one third of the range 
of ip in Fig.[|)]since the rest of the wave function can be ob- 
tained by symmetry: $(ip+4ir/3) = $(^+271-/3) = <&((/?). 
Further, we show the region that includes = — this 
point lies at ip = 7r/3 and 9 = tt/2. The hyperangular 
distribution tells us about the geometry of the system. 
At 6 = 0, for instance, the atoms form an equilateral 
triangle and at 9 — tt/2 they lie along a line. 

As R increases, a bound channel function collapses like 
i?" 1 to the region around = 0, as shown in the first 
row of Fig. [HI for a — 500 a.u.. The channel function oc- 
cupies a smaller and smaller region of the (ip, 9) plane as 
R increases, so must grow in amplitude to preserve nor- 
malization. The second row of Fig. [^represents the low- 
est continuum channel function for a — 500 a.u., which 
is the second channel for the system. Therefore, the 
channel function must approach the A = hypersphcri- 
cal harmonic (a constant function of the hyperangles) as 
R — * oo, and it must have one node (careful inspection 
of the figure near 9 = shows that $ < 0). In Fig. 
[SI the harmonic character is easily seen since the chan- 
nel function is spread nearly uniformly over the entire 
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«=1000a.u. ff=5000a.u. fi=10000a.u. 




FIG. 6: = + channel functions &(R; fi) for fixed values of the hyperradius R. The first and second lines show the bound 
channel (BC) and the lowest continuum channel (CC) for a = 500 a.u.. The third line shows the Efimov channel (EC), 
demonstrating that it is insensitive to changes in R. 



hyperangular plane. It is also apparent that there is a 
substantial amplitude at the two-body coalescence point 
and that this amplitude shrinks with growing hyperra- 
dius just as for the bound channel. 

The last row of Fig. shows the channel function for 
\a\ — oo. These channel functions share characteristics of 
both the bound and continuum channels. That is, they 
have large amplitude at the two-body coalescence points, 
but they are also spread over the entire plane. These 
channel functions, however, have one important charac- 
teristic that the others do not have: they are nearly scale 
invariant. That is, they look essentially the same at ev- 
ery value of R. This scale invariance is, in fact, one of the 
defining properties of the Efimov effect. These channel 
functions also show why the Efimov effect is so tenu- 
ous — the three bodies must maintain a delicate balance 
between the two-body and three-body parts of configu- 
ration space over ever-growing volumes of space. If there 
are two-body bound states present, this Efimov channel 
function will aquire nodes closely resembling those of the 
continuum channel. 

In Fig. \7\ we compare the channel functions for the 
weakly bound channel for a — 500 a.u., a — 1000 a.u., 
and \a\ = oo at R — 10000 a.u.. This figure shows 



the channel function on a logarithmic scale (in the z- 
direction) to emphasize the derealization in the Efimov 
case (\a\ = oo). At this value of R, the derealization 
grows with a, while all cases have large amplitude at the 
two-body coalescence points. 



C. Nonadiabatic couplings 

The asymptotic behavior of the nonadiabatic coupling 
for finite a has been deduced analytically for all combi- 
nations of bound and continuum channels [Tdj . The cou- 
pling to the Efimov channel, however, was only briefly 
discussed [2lJ, and that was only for the diagonal cou- 
pling element Q vv . We have found that the nonadiabatic 
coupling is also modified due to Efimov physics in the 
range ro -C R <C \a\. In Table [I] we show the leading 
term in the asymptotic expansions of P vv i and Q vu > cou- 
plings deduced in Ref. |T|| which are expected to be valid 
for R ^S> \a\. Here, we will restrict our analysis to cou- 
plings between s-wave channels and the Efimov channel. 
The numerically calculated couplings were obtained for 
values of the interaction strength near D — Dji (Fig. ^| , 
where the system has one bound state for a < and 
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FIG. 7: Channel function for the weakly bound channel at 
R = 10000 a.u. for a = 500 a.u. and a = 1000 a.u.. The 
|o| = oo channel function is also shown for comparision. Note 
that the vertical scale is logarithmic. 



two bound states for a > 0. In this case, therefore, we 
investigate the behavior of the nonadiabatic coupling in 
a situation similar to crossing an s-wave Feshbach reso- 
nance 0, EJ El 0- El , which involves the formation of a 
weakly bound molecule for a > 0. 
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BC(Z 


= 0) 


-BC(Z = 


= 0) 


R- 1 


R- 2 


BC(Z 


= 0)- 


-CC(Z 


= 0) 


R - 5 /2 


R- 1 ' 2 


CC(/ 


= 0) 


-CC(/ 


= 0) 


R- 2 


R- 4 



TABLE I: The leading term for the asymptotic expansion of 
the nonadiabatic coupling involving bound channels (BC) and 
continuum channels (CC). The diagonal coupling Q vv for BC 
is proportional to R~ 2 , while for CC it is proportional to R~ 4 
(see Ref. |H|). 

We plot in Figs. andllOlthe numerically calculated 
nonadiabatic couplings on a log- log scale in order to high- 
light the expected power law behavior. Figure [S] shows 
the couplings between an s-wave bound channel and the 
Ehmov channel for finite a (dashed lines) and \a\ = oo 
(solid line). For a > 0, and R 3> a, the couplings behave 
as couplings between bound channels, given in Table [I] 
For a < 0, they behave as couplings between bound and 
continuum channels. As \a\ — > oo, the couplings in both 
Figs.Efa) and (b) approach a limiting value given by the 
\a\ = oo result. In the range r <C R -C \a\, indicated in 
Fig. |H| by the vertical dot-dashed lines, the couplings are 
modified due to Efimov physics. The \a\ = oo coupling 
at the other poles of Fig. ^ will have the same behavior 
for i? 3> ro, but the short-range behavior will differ. 

Figure El shows the coupling between the Efimov chan- 
nel and the lowest continuum channel for a > and the 
next continuum channel for a < 0. The a > curves 
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FIG. 8: The nonadiabatic couplings between a bound channel 
and the Efimov channel for finite (dashed lines) and infinite 
(solid line) scattering lengths. For a > (a) the Efimov chan- 
nel is associated with the weakly bound channel; and for a < 
(b) with the lowest continuum channel. 



in Fig. EEa) show the characteristic peak around R ~ 3a 
predicted in [Tl| responsible for three-body recombina- 
tion at ultracold temperatures. For a < 0, there is one 
node at the dip in Fig. Elb) (since the absolute value was 
taken for this plot). As \a\ increases, the node moves to 
infinity, and the coupling for the Ehmov case emerges. 

The asymptotic behavior of the dia gona l correction in 
Fig.^Jalso confirms the results in Ref. [n| for finite a. In 
Fig. llUl thc diagonal coupling for \a\ = oo is proportional 
to i?~ 3 , which lies between the R~ 2 expected for a > 
and R~ A expected for a < (see Table 0. 

Our calculations confirm all of the finite \a\ asymptotic 
analysis of Ref. • Table [H] summarizes the cases for 
coupling with the Ehmov channel when \ a\ — oo. It shows 
the leading term in the expansion of the nonadiabatic 
couplings for three identical bosons and are expected to 
be valid in the range ro <C R <C \a\, as with the other 
results related to Efimov physics. Since the continuum 
channels are also modified due to the Efimov physics (see 
Fig. |3J, the couplings between bound channels and the 
continuum channels are expected to change. In fact, only 
the coupling between s-wave bound channels and contin- 
uum channels with Iq = are different from the finite a 
case. We have included this coupling in the last line of 
Table ITT1 for completeness. 
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FIG. 9: The nonadiabatic couplings between the Efimov chan- 
nel and a continuum channel for finite (dashed lines) and in- 
finite (solid line) scattering lengths. For a > (a) the Efimov 
channel is associated with the weakly bound channel; and for 
a < (b) with the lowest continuum channel. 



FIG. 10: The diagonal nonadiabatic coupling for the Efimov 
channel for finite (dashed lines) and infinite (solid line) scat- 
tering lengths. For a > (a) the Efimov channel is associated 
with the weakly bound channel; and for a < (b) with the 
lowest continuum channel. 
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EC 


-BC(Z = 0) 




R -3/2 
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-CC(Z = 0) 




R- 2 


R- 3 


BC(Z = 


0)-CC(Z = 


0) 


R -3/2 


R -5/2 



TABLE II: The leading term for the asymptotic expansion 
of the nonadiabatic couplings for \a\ = oo. The diagonal 
coupling Q„„ for the Efimov channel is proportional to _R -3 . 



Comparing Tables [I] and [HI we see that the Efimov 
channel shares characteritics of bound and continuum 
channels. The coupling between the Efimov channel and 
bound channels (i?~ 3 / 2 ) falls faster than the coupling 
between bound channels but slower than the cou- 

pling between bound and continuum channels (i?~ 5 / 2 ). 
The coupling between the Efimov channel and contin- 
uum channels {R~ 2 ) falls faster than the coupling be- 
tween bound and continuum channels (i?~~ 5 / 2 ), but, is 
the same as the coupling between continuum channels 

(ir 2 ). 

From Figs. 00 E3 and ^3 it is clear that the couplings 
for finite a also exhibits universal aspects. This universal 
behavior would better seen by plotting |a|P„„' (R/\a\) and 
\ a \ 2 Quv'(R/\ a \)- In such a plot the couplings for different 
values of a look all the same for R/\a\ 3> 1. 



IV. SUMMARY 

We have investigated several fundamental aspects of 
Efimov physics in the hyperspherical adiabatic represen- 
tention of systems of three identical bosons. We have 
calculated effective hyperspherical potentials for very 
large to infinite values of the two-body s-wave scattering 
length. These calculations were obtained in a regime not 
accessible in previous calculations, allowing us to test as- 
sumptions related to the Efimov effect, namely, the emer- 
gence of the long-range three-body effective interactions 
in the range where ro -C R <C \a\. Our results show that 
despite the more restrictive limits discussed here for the 
effective potentials, the assumption that the long-range 
attractive potential holds in the range ro ^ R ^ \a\ gives 
a good qualitative description of the Efimov physics. In 



fact, for 



> 



ro, the unrestricted result for the number 



of Efimov states seems to give a better description, even 
though the effective potential does not resemble the Efi- 
mov potential. Additionally, we have shown that as we 
increase the number of two-body bound states the effec- 
tive potential seems to converge to a limiting curve. 

We have also obtained the behavior of many other 
quantities affected by Efimov physics: the continuum ef- 
fective potentials, all of the nonadiabatic couplings with 
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the Efimov channel, and the channel functions them- Acknowledgments 
selves. As expected, the channe functions are essentially 
unchanged with R. Moreover, the nonadiabatic couplings 
reveal that the Efimov channel behaves in some ways like 
a bound state and in others like a continuum state. We 
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